Hydrodynamics of the quantum Hall smectics 
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We propose a dynamical theory of the stripe phase arising in a two-dimensional electron liquid 
near half-integral fillings of high Landau levels. The system is modelled as a novel type of a smectic 
liquid crystal with the Lorentz force dominated dynamics. We calculate the structure factor, the 
dispersion relation of the collective modes, and their intrinsic attenuation rate. We show that thermal 
fluctuations cause a strong power-law renormalization of the elastic and dissipative parameters 
familiar from the conventional smectics but with different dynamical scaling exponents. 
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Planar arrays of interacting lines or stripes have be- 
come a paradigm for many different physical systems, 
including domain walls in magnets, layered supercon- 
ductors, biophysical systems, liquid crystals, and charge 
density waves (CDW) §. A unidirectional CDW or the 
stripe phase was also predicted || to form in GaAs two- 
dimensional (2D) electron systems when the occupation 
of the third or a higher Landau level is close to i || . This 
prediction was attested by a recent discovery of a dra- 
matic magnetotransport anisotropy Q in the indicated 
range of filling factors. The easy (low resistance) and 
hard (high resistance) current directions are thought to 
be along and across the stripes, respectively. In order to 
have a better foundation for the magnetotransport stud- 
ies, one has to first understand the intrinsic dynamics of 
the stripe phase, without intervening disorder effects ||. 
Recent investigations |^-^] revolved around the instabil- 
ity of the stripes against the 2fcF-modulation along the 
direction of the translational order, which would trans- 
form the system into a highly anisotropic Wigner crystal. 
According to MacDonald and Fisher [|| , such an instabil- 
ity appears only at temperatures T below 1 mK. Here we 
would like to address another interesting regime of rel- 
atively high T where the 2kp modulation is absent and 
moreover, the sole periodic modulation is due to the main 
CDW harmonic with wave vector q$. In classical terms, 
the electron density at the topmost Landau level is of the 
form n(r,t) + Re#(r,i), where * = |$|e i *°0 ! - ,u ) is the 
CDW order parameter. (We chose the x-direction to be 
perpendicular to the stripes). Although oversimplified, 
this classical picture correctly identifies the low-frequency 
long-wavelength degrees of freedom: the coarse-grained 
component n of the total density, and the CDW phase 
(or the displacement field) u. It is natural to expect that 
at sufficiently small w and q the dynamics of u and n is 
governed by a certain hydrodynamic theory. The formu- 
lation of such a theory is the subject of this Letter. 

We begin by writing down an effective Hamiltonian, 
from which we derive the hydrodynamic equations of mo- 



tion, the spectrum of the collective modes, and the struc- 
ture factor. Next, we study the renormalization of the 
bare parameters of the theory by thermal fluctuations. 
Our results can be verified by measuring the transmis- 
sion of the micro- and surface acoustic waves as well as 
inelastic light scattering techniques provided the samples 
are of high enough mobility so that the disorder effects, 
which we ignore here, are not important. 

Effective Hamiltonian. — Our first step is to construct 
an effective Hamiltonian for n and u, which is (a) rota- 
tionally and translationally invariant and (b) reflects the 
specific properties of the system, such as the long-range 
Coulomb interaction. Very useful in this process is the 
similarity Q to smectic liquid crystals Q] . The resultant 
form of the effective Hamiltonian is as follows: 

H=^J al 2 r[YE 2 (u) + E(u)A + K(V 2 u) 2 

+ X~ 1 Sn 2 + <5>n + 2CE(u)6p + pv 2 }, (1) 

where E(u) — V x « — |(Vji) 2 in the usual rotationally 
invariant strain ||, Y and K are the compression and 
bending elastic moduli, x 1S t ne compressibility []lo|| , A is 
an auxiliary parameter (counter-term), needed to guar- 
antee the condition (V^u) = 0, p = mn is the mass 
density, Sp = p — po is its deviation from the equilibrium 
value po — mriQ. <f>(r) = J d 2 r'Sn(r')U(r — r'), is the 
electrostatic (Hartree) potential, with U(r) = e 2 / nr at 
large r (see more details in Q). The penultimate term in 
Eq. ([[]) accounts for the dependence of the CDW period- 
icity on no, with C = (Y / m)d In qo / dn$. Finally, the last 
term is the correction to the kinetic energy, with v being 
the velocity of the electron fluid. This term is mainly a 
bookkeeping device: it vanishes after the projection on 
a single (topmost) Landau level but enables us to derive 
the equations of motion via the standard Poisson-bracket 
method 0H- 

Equations of motion. — The hydrodynamic fluctua- 
tions of our system are governed by the equations 
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F.j = + VjipvkVj) - ptu c e jk v k + VjW— 

+ L jk (iV)v k + pV fe (x^Sp + <S> + CE) = Q, (3) 

and the continuity equation p + V(pv) = 0. The no- 
tations here are as follows: W = x — u, Ej k is the unit 
antisymmetric tensor, and lo c — eB/mc is the cyclotron 
frequency. Functions Lj k (q) describe dissipation. As in 
3D case ||, they can be parametrized by viscosities rji, 



+ Sjx$kx[mq 2 + im - 4% + m - 2^)q x \ 

+ ("3 -Vi + V5)qx(q 3 S kx + q k 5 jx ). 



(4) 



(There is no analog of 7/2 in 2D) . One more dissipative co- 
efficient, A p , [Eq. (||)] describes the permeation, i.e., the 
mass transport across the stripes. At relatively high T, 
such that |*| < n , X p ~ (n 3 q^/ p Q uj 2 )(n /\^\) 2 . Finally, 
£ a 's in Eqs. (||) and (||) are the Langevin noises, which 
satisfy the fluctuation-dissipation theorem, 

(C(l)Cfc(2)> = 0, l(2) = {r 1(2)j t 1(2) }, (5) 
(C„(1)C™(2)> - 2k B TX pPo S(l - 2), (6) 
<0(1)C*(2)> - 2k B TL jk (iV)S(l - 2). (7) 

The main difference of our model from the conven- 
tional smectics is the presence of the strong magnetic 
field. As we will see below, it drastically changes both 
the linearized and nonlinear dynamics of the system. 

Field-theory description. — Our ultimate goal is to cal- 
culate various correlation functions, e.g., the dynamical 
structure factor, 



where No is the total number of electrons at the topmost 
Landau level, and (...) stands for the thermal averaging, 
or equivalently, the averaging over the Langevin noise. 
The latter can be done in a systematic way by means 
of the Martin-Siggia-Rose (MSR) formalism p3| , whose 
field-theoretic version was previously applied to the con- 
ventional smectics by Kats and Lebedev [[l2| . 

The implementation of the MSR method begins with 
enforcing the equations of motion by 5-function-type 
weight factors, S(F a — ( a ), in the path integral over F a . 
These 5-functions are then represented by integrals of 
exp[ip a {F a — £ a )] over auxiliary dynamical variables p a , 
which enables one to average over the Gaussian noise Co- 
Finally, one changes the integration variables from F a 
to Vi = pvi/po and u. If the denote by F™ 1 the non- 
liner terms in F a , and introduce vectors (jr = {Vt,pt}, 
yt = {u, V x , V y }, = {p u ,Px,Py}, then the resultant 
action becomes 



A = -i j dt J d 2 r Pa F" l {u, v) - In J + A , (9) 



An = - 



where J = | det dF a /dVb\ is the Jacobian and Go is the 
(bare) propagator of the following block-matrix form 
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(11) 



Two particular components of Go , G uu and G 
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G uPx (q, w) 



G uu (q, u) 



Q{q)s 2 — iu;D(q) — lo 2 - 



(12) 



k B T 
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[G uP:c (q,ui) -G MPx (q, -w)], (13) 



will play an important role in the later discussion. Here 
we introduced the notations Q(q) = (Yq 2 + Kq 4 )/p , 
w p(<?) = n q 2 [U{q) + X _1 ]/ m > h = %<? 2 /Po, c = q x /q = 
cos 9, s — q y /q~ sin#, and 



^(q) = ^3 + c 2 s 2 (i>i - Av 3 + z> 4 - 2i> 5 ). 



(14) 



Equations (|12j) and (|13|) are obtained by neglecting u>, X p , 
and C compared to quantities proportional to "large" 
frequencies u> c and u> p . 

Harmonic theory. — If only the quadratic part Aq of 
the full action is retained, then the correlators of the ve- 
locity fields V a are given by the components Gv a v b of Go 
and can be found after straightforward albeit tedious al- 
gebra. The continuity equation enables us to relate these 
correlators to the structure factor. This way we get 



S(q, W ) = M^ / dte M *(p q (0p- q (0)}, (8) 5(q,a;) 
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k B Tq 2 (Qs 2 - iujv) 



mui Qlo 2 (s 2 + \ p v) — iujuj 2 Oi — lo 2 lo 2 ' ^ ^ 



for lo <C lo c and, in particular, 

q 2 U! 2 X p + c 2 v 3 + s 2 i> 4 



5(q,0) = 2fc B T- 



s 2 + X z 



(16) 



The poles of S(q,oj) correspond to the collective 
modes. Two of them (magnetophonons) are gapless: 



w m (q) ~ sin 20 
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u p {q)q 

2UJr 



*(q) (17) 



and its counterpart u) = — uj* n (c\). Their attenuation rate 
a(q) is given by 

a(q) = A P Q + oj- 2 [uj 2 p (q)D + Q(D 3 c 2 + D 4 s 2 )}, (18) 

up to terms proportional to C and higher powers of Vi and 
Xp. The other two collective modes (magnetoplasmons) 
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have a large gap oj c . Thus, the number of the hydrody- 
namic (gapless) modes coincides with the number of the 
hydrodynamic variables (n and u) as it should p3 . 

For Coulomb interaction uj p (q) oc q 1 ^ 2 ] hence, within 
the harmonic theory the magnetophonons have oj m (q) oc 
q z l 2 dispersion, similar to that of the Wigner crystal 
but with the ^-dependent prefactor. Propagating mag- 
netophonon modes exist as long as 9 is not too small so 
that Rew m (q) » a(q) in Eq. (|l7|); otherwise, they are 
replaced by the two overdamped modes: 



wfast(q) ^ -ia(q), w s iow(q) 



a(q) 



q x q y m P \ Y_ 
q u c ) Pa' 



Below we will see that at small enough q the magne- 
tophonon dispersion relation and their damping are sig- 
nificantly modified by the anharmonisms. 

Renormalization. Our next step is to calculate the 
propagator G = (Gq 1 + of the full theory, Eq. (§), 
treating the previously ignored nonlinear terms as per- 
turbations. We expect that G has the same form as 
Go but with Y, K, and other parameters replaced 
by oj and g-dependent (renormalized) values Y R (q, w), 
K R (q,uj), etc. The perturbative corrections to Y, K, 
and C are determined by the self-energy component 
E uPa .(q, lu). Similarly, the corrections to Vi's are deter- 
mined by Sp.jp,. (q, <jj). To the lowest order in T they are 
given by the diagrams shown in Fig. |l|. 




FIG. 1. Self-energy diagrams. 

Doing the power counting with the help of Eqs. (|1£ 
(|l|), and the relation k x ~ ^{K/Y) 1 / 2 , where k is the 
loop momentum, we quickly discover that S„ Px and S Pj . p( . 
are infrared-divergent. To resolve this problem, we uti- 
lize the renormalization group (RG) procedure formu- 
lated in d = 3 — e spatial dimensions [one x-dimension 
plus (d — 1) "_L" -dimensions]. Each step of our RG trans- 
formation consists of integrating out the cylindrical shell 
Ae~ l < k± < A in the loop diagrams of Fig. |[ followed 
by rescaling of the momenta k = k'e~ l to restore the ul- 
traviolet cutoff A ~ qo on k± . It is convenient to rescale 
the u-field as well it(r) = u'(r')e l to preserve the struc- 
ture of E(u). We did not find it necessary to rescale the 
frequencies or impose a cutoff on k x The one- loop 

RG equations are given by 



d+2 



K. 



7l 9= (3-% 
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dl 



2dv{ 
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2(d-l) 

,2 



K ( l ) c (d+2)l + Yl 

cj2(Ae-') po 4 



(19) 
(20) 
(21) 

(22) 



where g = ksTY 1 / 2 K^/ 2 / S d -i is the dimension- 
less coupling constant, with Sd-i = (47r)~s~r (^rp) /2. 
Equations for other viscosities are similar to (||^) and are 
not shown. 

The solutions of the RG equations (|l9|-|2^) are as fol- 
lows: 



g(l) = g* + 



9o - g* 



D = 1 



50 [p (3-d)Z 



D g„ ■ 

Y(l) = Y e dl D-~< Y , C{1) = C Q e dl D^ Y , 



1], 4{l) = e 
„d-l 



2,11 



"1(0) 



2(d-l) W a(Ae-') Po 



DXd+2) . 



1Y 



J d+2' 



K(l) = K eS d -^ l D^ K , 1K = 



d + 2' 



(23) 



(24) 



As one can see, for small e = 3 — d, g flows to a weak 
coupling fixed point g* = 16(3 — d)(d — l)/(rf + 2), jus- 
tifying our one-loop RG. The renormalized values of the 
parameters of the harmonic theory are found by rescaling 
back to the original coordinates, Y R (ui, q) = Y(l)e~ dl , 
K R (u,q) = AT(/)e-( d - 2 )', is R (uj,q) = v 3 (l)e- dl , etc., 
where I is the smallest of the three cutoffs l x , l±, Z w to 
be found from equations K(l x )K A = q 2 e 2l:c Y(l x ), A = 
q±e l± , and uj 2 c uj = u 2 (Ae- l -)p 3 (L)A 2 e- {d+ ^ L . These 



I 

equations give rise to the characteristic crossover length 
and frequency scales £ x = (K^ /F ) 1/2 / '(k B T) 2 , £ y = 
(Kl/Yo^/ksT, and u y = (K Q / Po )V 2 w p (^)/lo c $. 

Taking Eq. (p4| ) for the face value, we obtain jy = 
7x = \ in 2D (e = 1) and discover the following three 
types of asymptotic behavior: 

Y R ~y (^) 1/2 , K R ~Koti v q v r 1/2 , 

V R K o V R _ a * Y o t 2 <\ 

3 Po^egtovfv) 3 / 4 ' 1 P^ y {q y (i y )y± 1 j 
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for q x « S-^qyty) 3 ' 2 , q y < ^ , U « U y (q v Z v p 4 , 



Y 1 



(26) 



3 P0^Cy(^6) 1/2 ' * P0^ y ((? 2; ^) 7/6 

for efe < f" 1 , <7y < C^H^^x) 273 , w < Wj,(gx^) 3/2 , and 

(27) 



Vt 2/3 ^ 1/3 , 





2 /9, .7/9 



for w > c^((7 y £ y ) 9 / 4 , w > Wj,(g x ^ x ) 3 / 2 . Here ai = 
(1 — r) 2 , a4 = r 2 , 05 = (1 — r)r, and r = pqCq/Yq. Thus, 
0,104 = a 2 , which entails the relation jnj A7/1A774 = 
(A775) 2 derived earlier for the conventional smectics. 

Discussion. Our one-loop-level results for Y R and K R 
turn out to be exact jL7| for the static limit [Eqs. ( p^ ) 
and (p6|)]. It is possible that this success transcends to 
the dynamics, in which case the scaling exponents for 
viscosities in Eqs. (p5|-p7j) would be exact as well. 

Within the domain of anomalous hydrodynamics sum- 
marized by Eqs. (|25|-[2^), the magnetophonon dispersion 
relation becomes 



Re Wm (q) ~ sc V^ x q)^^p- J*L, 



Imw m (q) ~ S 2 c 5 / 6 (Cx?) 11/6 (l-2r) 




(28) 
(29) 



for a finite fixed 9 and q <C £ x 1 - In contrast to the case 
of conventional smectics the scaling of Imw m (q) 

is cut off by l x [Eq. @] not l w [Eq. @]. 

Strictly speaking, the smectic behavior in 2D can per- 
sist only up to the length scale smaller than the aver- 
age distance £<j ~ exp(£' £ ;/2fc_BT) between dislocations, 
Ed ~ i^o being the dislocation energy |Q. The anoma- 
lous hydrodynamics can be observed if £ x < £,d, which is 
satisfied at low T where £d is exponentially large Jl9| . 

In conclusion, we formulated a novel long-wavelength 
low-frequency effective theory of the stripe phase aris- 
ing in a two-dimensional electron liquid near half-integral 
fillings of high Landau levels. Our theory applies at rel- 
atively high temperatures and in the clean limit. We 
demonstrated that the collective mode properties of the 
system exhibit nontrivial power-law scaling, which can 
be verified by microwave or surface acoustic wave mea- 
surements at a finite wave vector or by inelastic light 
scattering experiments. Our future plans include apply- 
ing the proposed hydrodynamic approach to the problem 
of the magnetotransport. 
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